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Method and apparatus for determining hysteresis 

FIELD OF THE INVENTION 

[0001] The invention relates to process automation systems. 

BACKGROUND OF THE INVENTION 
5 [0002] In a process automation system, a control loop typically con- 

sists of a process, a measurement, a controller and an actual control element, 
such as a valve, and a related device, such as a valve controller (positioner) 
and an actuator. An optimal process control depends on how appropriately all 
these components function. 
^0 [0003] In the processing industry, such as the pulp and paper indus- 

try, and petroleum refining, petrochemical and chemical industries, various 
control valves mounted in the plant pipe system control material flows in the 
process. A material flow may contain any fluid material, such as flowing sub- 
stances, liquors, fluids, gases or vapours. At its simplest, the control valve may 
15 be a manually-controlled mechanical valve. Usually, the valve is provided with 
a valve controller and an actuator. The valve controller and actuator adjust the 
position of a control valve according to the control input (e.g. pneumatic or 
electric control input) received from the process control system. 

[0004] Figure 1 exemplifies a functional block diagram of a control 
20 valve. A valve controller (i.e. positioner) 10 controls the travel/position (h) of 
the valve by means of the torque generated by an actuator 11. Position infor- 
mation (h) is provided as feedback from the actuator 11 or valve 12 to an ad- 
der arranged at the input of the valve controller. The function of the valve con- 
troller is mainly based on an error (e) between an input signal u (control signal 
25 from the process) and a feedback position (h). The valve controller 10 mini- 
mizes this error by a control algorithm, such as a state or PID algorithm. This 
control algorithm is tailored for each valve and, rf necessary, it can be tuned 
during the Installation or operation. The tuning may include changing gain pa- 
rameters. It is also feasible to use one or more additional feedbacks in the 
30 valve controller 10, such as speed or pressure feedback from the actuator's 
cylinder, to obtain a more balanced and accurate control function of the valve 
position. 

[0005] Figure 2 illustrates a typical model of a process control loop 
which controls one control valve 22 and therethrough one material flow in a 
35 process. The control valve 22 may be similar to the one shown in Figure 1 . for 
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instance. The process control loop includes a process controller 21 . which is 
provided with a control algorithm, which produces a control signal u for control- 
ling the control valve 22 according to a set point r (which is obtained from a 
control room computer, for example) and a feedback process variable y. The 

5 control algorithm may be any algorithm that is used in control systems, such as 
PID, PI or P control. The control signal u fed to the control valve 22 controls 
the valve position and travel and thus the material flow in the process. A de- 
sired process variable y is measured by a measurement transmitter 24 and it is 
compared (block 20) to the set point r of the same process variable to produce 

1 0 an error signal e2, which is fed into the process controller 21 . The process con- 
troller 21 changes the control signal to minimize the control error 2e. The proc- 
ess control error typically results from changes in a positioning error and proc- 
ess interference. 

[0006] The valve and its auxiliary devices often constitute the weak 
15 link in the control loop since they are the only moving parts. This movement 
causes problems, which decrease the capacity of the control loop. To avoid a 
backlash resulting from mechanical adjustments, the valve, actuator and valve 
controller/positioner have to be provided with mechanical tolerances that are 
sufficiently tight. As a result of the backlash, the valve movement does not fol- 
20 low the control signal accurately but deviates from it. The influence of the 
backlash becomes apparent in particular when the valve control direction and 
thus the valve's direction of movement are reversed. In that case, the control 
signal value keeps changing for a while until the measured output signal starts 
to change noticeably. This is also known as the dead band of control. In addi- 
25 tion to the backlash in an actuator or positioner, this phenomenon may result 
from sticking of the valve or other mechanical factors, such as initial friction. 
The backlash between mechanical parts naturally increases as the parts wear. 

[0007] The backlash and other error factors cause hysteresis be- 
tween the control of the process device, such as a valve and/or its auxiliary 
30 devices, and the measured response. This is illustrated in Figure 3. Straight 
line 31 illustrates an ideal relation, i.e. characteristic curve, between the control 
u and the measurement (output) y, such as valve position. The real depend- 
ency between the measurement and the control is illustrated by characteristic 
curve 32. As appears from Figure 3, due to the backlash and any other factors, 
35 the upward control (increasing u) has a characteristic curve 32A different from 
that of the downward control (decreasing u), which has characteristic curve 
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32B. The difference between the curves represents hysteresis in the control of 

the process device. 

[0008] In some cases, the controllers are provided with automatic 
bacl<lash compensation, which attempts to take the mechanical non-ideality of 
5 the device into account always when the control direction is reversed. This ap- 
proach is described in US patent 5,742,144, for example. Approach of this kind 
is good in theory but limited in practice since the backlash and hysteresis vary 

due to different factors. 

[0009] The information on hysteresis and backlash is, however, im- 
1 0 portant to the tuning of the control circuit. It also gives useful information on the 
condition of the process device, such as a valve and/or its actuator or valve 
positioner. If hysteresis or backlash increases significantly, service measures 
can be taken to fix the matter. 

[0010] A typical way of detecting the hysteresis or backlash of an 
15 actuator Is to switch the controller to a manual controlling mode and perform a 
sequence of step tests. In that case, the actuator is driven to the same position 
from different directions, in which case any differences between the control 
and the response due to backlash or hysteresis are found out by means of 
measurements. Another typical way is to drive the actuator back and forth over 
20 the whole control area and to estimate backlash and hysteresis from the 
measurement results. In the case of a valve actuator, for example, the valve is 
driven from the closed position to the open position and back to the closed po- 
sition. A problem associated with these solutions is, however, that they are 
separate tests that need to be carried out when the process is Interrupted or 
25 the process device to be examined is bypassed or detached from the process. 
Similar tests that are performed on the valve positioner are described In lEC 
(International Electrotechnical Commission) standard 61514, Industrial process 
control systems: Methods of evaluating the performance of valve positioners 
with pneumatic outputs, first edition, 2000-04. 
3Q [0011] WO 01/11436 discloses a method and an apparatus which 

statistically detennine estimates for one or more process control loop parame- 
ters for the device or the control loop that is active in the process control envi- 
ronment. Such parameters include friction, dead band, dead time, vibration or 
backlash. In the method, one or more signals are always measured in the 
35 process control loop when the process control loop is connected to the on-line 
process control environment. The measured signal is stored as signal data, 
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after which a number of statistical analyses are carried out on the stored data 
to determine the desired parameter estimate. An advantage of this solution Is 
that the process device does not need to be removed from process or the con- 
trol loop bypassed for the test. 
5 [0012] In practice, the on-line detemfiination of hysteresis or back- 

lash is sensitive to process interference as well as inaccurate. Furthermore, it 
usually requires statistical calculation methods, matrix calculation, mathemati- 
cal functions, etc. 

BRIEF DESCRIPTION OF THE INVENTION 
-10 [0013] The object of the invention is to provide a new method and 

apparatus for determining hysteresis or backlash of a process device in normal 
run in a process in a relatively simple and accurate manner. 

[0014] The object of the Invention is achieved by the method and 
system described in the attached independent claims. Prefenred embodiments 
15 of the invention are described in the dependent claims. The present invention 
comprises collecting sample pairs of a signal representing the control of a 
process device in normal run and a signal representing its output. The signal 
representing the output can be, for example, a process variable measured in 
the process. It may also be the set point of a process variable. Each collected 
20 pair of control/output samples comprises an average control input and a 
measured average output, which have been calculated for a predetermined 
collection period. In an embodiment of the invention, a momentary sample pair 
is typically taken from the signal representing the control input and from the 
measured signal representing the output at certain intervals, the interval being 
25 preferably in the order of one or more seconds. Such data is often called 'sec- 
ond-level data' or 'seconds data'. As in one embodiment of the invention, an 
average sample pair is then calculated from the momentary sample pairs (e.g. 
seconds data) taken during the collection period. The collection period is pref- 
erably in the order of one or more minutes, and thus the mean values are often 
30 called "minute mean values'. 

[0015] Before calculating characteristic curves, sample pairs suit- 
able for further processing are screened from the collected raw data on the 
basis of the magnitude and direction of the relative change of the average con- 
trol input. This screening is of great importance to the accuracy of the method. 
35 In normal run, the control circuit continuously controls the process device so 
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that the control zigzags in both directions. In addition, various kinds of interfer- 
ence occur. By screening suitable sample pairs for further processing on the 
basis of the magnitude and direction of the relative change of their control in- 
put, corner or turning points of the control situation where the control direction 
5 changes can be found out. As the control direction changes, the measured 
output may fall in a vague area due to hysteresis or backlash and thus cause 
errors in the final calculation of characteristic curves. The screening according 
to the Invention ensures that the samples selected for calculation represent a 
stable situation either in upward control or in downward control. In an embodi- 
10 ment of the invention, this selection or screening is canried out by calculating 
the change in the average control input in relation to the average control input 
of the previous sample pair and the direction of change for each sample pair. If 
the calculated change is smaller and has occurred in the same direction as the 
change calculated for the control input of the previous sample pair, the sample 
1 5 pair is selected for further processing; otherwise the sample pair is rejected. 

[00161 According to the invention, the sample pairs are also 

grouped into a first group and a second group according to the direction or 
sign of the relative change of the average control Input of each sample pair. 
These two groups represent upward control (increasing control value, positive 
20 sign) and downward control (decreasing control value, negative sign). The 
grouping can also take place before, during or after screening. By means of 
the screened sample pairs of the first and the second group, a first and a sec- 
ond control/output characteristic curve representing ascending and descending 
control are calculated. Finally, the hysteresis or backlash of the process device 
25 is determined on the basis of these characteristic curves as a distance be- 
tween the first and the second characteristic curve at least at one point. In 
practice, the calculation of a characteristic curve comprises calculation of at 
least a few points for each characteristic curve. The detennination of hystere- 
sis or backlash means that the distance parallel with the control axis Is deter- 
30 mined at least at one point detemiined this way. 

[00171 It is important for the operators of an Industrial process to 
know how reliable the hysteresis or backlash obtained in the manner according 
to the invention is. The reliability varies with different process conditions. In an 
embodiment of the invention, a reliability value is also calculated for the deter- 
35 mined hysteresis or backlash by a function, which includes the following infor- 
mation as parameters 
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the number Ny1 of average sample pairs belonging to the first, as- 
cending characteristic curve, 

the number Ny2 of average sample pairs belonging to the second, 

descending characteristic curve, 
5 the number Noy1 of sample pairs which belong to the first charac- 

teristic curve but are below the second characteristic curve, and 

the number Noy2 of sample pairs which belong to the second char- 
acteristic curve but are above the first characteristic curve. 

[0018] This determination of reliability is based on the fact that the 
10 higher the number of the measurement pairs which are used for calculating a 
certain characteristic curve and are below the second curve, the more unreli- 
able the result is. According to an embodiment of the invention, the reliability 
value is calculated by function L=max(0,1-Noy1/Ny1-Noy2/Ny2). in which case 
L may obtain values from 0 to 1, where L=1 is completely reliable and L=0 is 
15 completely unreliable. 

BRIEF DESCRIPTION OF DRAWINGS 

[0019] In the following, the invention will be described by means of 
exemplary embodiments, with reference to the accompanying drawings, in 
which 

20 Figure 1 is a functional block diagram of a control valve, 

Figure 2 illustrates a typical model of a process control loop, and 
Figure 3 is a graph illustrating a characteristic curve and hysteresis 

of the control loop. 

Figure 4 is a schematic functional block diagram illustrating an algo- 
25 rithm according to an embodiment of the invention for detemninlng hysteresis. 

Figure 5 is a data flow chart illustrating the algorithm shown in de- 
tection block 40 of Figure 4, 

Figure 6A illustrates minute mean values of measurement y. 

Figure 6B illustrates minute means values of control u as a function 

30 of time. 

Figure 7A illustrates characteristic curves fomned without the selec- 
tion according to the invention, 

Figure 78 illustrated characteristic curves formed using the selection 

according to the invention. 
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Figures 8 and 9 illustrate calculation of hysteresis from the charac- 
teristic curves, and 

Figure 10 is a graphical presentation where the upper image illus- 
trates control measurement pairs and a calculated characteristic curve, and the 
5 lower image illustrates a frequency curve. 

DETAILED DESCRIPTION OF THE INVENTION 

[0020] The present invention is applicable to all industrial processes 
and the like. The invention is not limited to any particular process device but it 
can be applied in various devices that control a process, such as valves, 
10 pumps, fans, heat exchangers, etc. 

[0021] The exemplary embodiments of the invention will be de- 
scribed using control valves and their auxiliary devices, i.e. valve controllers 
(positioners) and actuators, as examples. 

[0022] Figure 4 is a schematic functional block diagram illustrating 
15 an algorithm according to an embodiment of the invention for detemiining hys- 
teresis. The controller 21 is similar to the PID process controller shown in Fig- 
ure 2, for instance. It receives as input a set point r of certain type and a 
measured process variable or another measurement y, which represents the 
control response or output of the control circuit or process device to be moni- 
20 tored. The controller 21 provides control u, which is fed into the process de- 
vice, such as a valve controller, an actuator or a positioner. Block 40 repre- 
sents hysteresis identification according to an embodiment of the invention. In 
the example of Figure 4, block 40 collects samples from the control signal u 
and feedback measurement result y. The measurement result may be the 
25 measured value of a controllable process variable, such as the measured flow 
as valve output. Instead of the actual controllable process variable, the meas- 
urement result y may be any measured process variable that represents the 
output and correlates appropriately with the controllable process variable or 
output. Altematively, the measurement result may be replaced with a sample 
30 representing the set point r. This may be done particulariy In cases where the 
measurement value y has been noted to follow well the set point r. This is the 
case In the example illustrated in Figure 6A, for instance. 

[0023] In an embodiment of the invention, the hysteresis detection 
block 40 follows the algorithm illustrated by the data flow chart of Figure 5. 
35 Block 40 collects second-level sample data ys and Us on the process from the 
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control and measurement signals u and y of the control circuit. In the second- 
level data, a typical sampling rate is one second or a few seconds, but the in- 
vention is not limited to any particular sampling rate. This sample data is proc- 
essed in different stages and as a result, the calculation provides an estimate 
5 for the hysteresis and reliability of the process device or control circuit. The 
method can be divided into an on-line phase and an off-line phase, but this 
division is not necessary or it can be performed differently from what is de- 
scribed in this example. 

[0024] In the on-line phase, second-level measurements are col- 
10 lected from the process control circuits or from the controller 21 control u and 
measurement y using another collecting rate, as illustrated in phase 50 in Fig- 
ure 5. IVlinute mean values ym and Um are calculated from the second-level 
measurement values in step 51 and these are stored in a database in step 52. 
The minute mean value means that the mean values are calculated from sec- 
15 ond-level values collected during a collection period whose duration is ap- 
proximately one minute or in the order of minutes. The collection period over 
which the mean value is calculated may vary depending on the application. 
The controller speed or integration time Ti is typically taken into account in the 
selection of a mean value period so that the updating period Tp (collecting 
20 period) of mean values fulfils the condition Ti < Tp <5 x Ti. If the controller's 
integration time T| were short in relation to the updating period Tp. the controller 
21 would have enough time to perform several controls during one period Tp. 
To give an example, Ti is typically in the order of 10 to 20 seconds in flow 
control. [0025] If the automation system already includes on-line collection 
25 for another purpose, by means of which minute mean values, for example, are 
collected from the control and measurement signals of the control circuit, this 
information can be utilized as such for the purpose of the present invention. In 
that case, the existing operations would replace steps 50, 51 and 52 in Figure 
5. 

30 [0026] The minute-level (u, y) pairs stored in the on-line phase are 

retrieved later, for example once a day, once a week or at another interval, 
from the database for hysteresis calculation. The pairs suitable for hysteresis 
calculation are selected from the stored minute-level sample pairs Um and ym 
according to a certain procedure (step 53). The pairs are also divided into two 

35 groups. Unsuitable pairs are rejected (step 57). Two characteristic curves (54) 
are calculated from the selected pairs (Um, ym) for calculating hysteresis (step 
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55). The calculation 55 also includes a routine which estimates the reliability of 
the identified hysteresis. 

[0027] In the following, implementation of different phases of the al- 
gorithm according to the invention will be described in greater detail by exam- 
5 pies. 

Calculation of minute mean value 

[0028] In an embodiment, the second-level sample pairs (or sample 
pairs taken at another collecting interval) from the measurement and control 
signals y and u are collected into buffers ys(i) and Uso), where i = 1...N (the 
10 buffer length is N = 60 when the second/minute calculation is used). The min- 
ute-level mean values ym and Um or mean values with another storing fre- 
quency are calculated on the basis of the seconds data included in the buffers. 

Selection and division of pairs of control and measurement 

[0029] In an embodiment before calculating a characteristic curve, 

1 5 the pairs of control and measurement samples suitable for hysteresis calcula- 
tion are selected from the minute-level raw data. The selected sample pairs 
are further divided into two groups G1 and G2 depending on whether the con- 
trol circuit control goes upwards (the value of control u increases) or down- 
wards (the value of control u decreases). 

20 [0030] A pair (UmO). ymO) of minute mean values is selected for the 

calculation of characteristic curve if the following two conditions are fulfilled: 

K,(y)|<|AM„,(7-i)l 

sign{Au„ U)) = sign{Au„ (J - 1)) 

25 where 

UmG) is the average control input of sample pair j, 
YmG) is the average measured output of sample pair j, 
j is an integer index, 

AUmG)= UmG)-UmG-1) 

30 sign(AUm) calculates the sign of AUm- • 

[0031] In other words, a pair of measurement and control samples 
is accepted for hysteresis calculation only if the change AUmG) the minute 
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mean value UmG) of the control is smaller than the change calculated for the 
control of the previous sample pair and the change has occurred in the same 
direction as the previous change. 

[0032] The pairs of control and measurement samples are further 
5 divided into two groups on the basis of the sign of the control change AUm as 
follows: a sample pair is taken in group G1 when the sign of the change is 
positive (the control goes upwards) and a sample pair is taken in group G2 
when the sign of the change is negative (the control goes downwards). 

[0033] The selection and grouping carried out according to the prin- 
10 ciples of an embodiment of the invention are very essential to the accuracy 
and reliability of hysteresis determination. This can be examined by means of 
an example related to flow control. Minute mean values are collected from 
measurement and control signals of a control circuit for approximately 6.5 
hours according to the invention. Figure 6A illustrates minute mean values of 
15 measurement y and Figure 6B illustrates minute mean values of control u as a 
function of time. It can be seen from Figures 6A and 6B that the control signal 
u zigzags; a phenomenon that may indicate a fault in the operation of an ac- 
tuator. If pairs of control and measurement samples are fonmed from the min- 
ute mean values shown in Figures 6A and 6B and characteristic curves are 
20 drawn without the selection according to the invention, we obtain the result 
shown in Figure 7A. It can be seen that the pairs of control and measurement 
samples fill the whole hysteresis area and no clear curves are discernible. 
Next, we will examine an embodiment according to the invention where the 
pairs of control and measurement samples suitable for hysteresis calculation 
25 are selected and the pairs are divided into two groups. The selected samples 
pairs are shown in Figure 78, where the sample pairs of group G1 (control 
moves upwards) are marked with symbol o and the sample pairs of group G2 
(control moves downwards) are marked with symbol x. Furthenmore, two char- 
acteristic curves -*- have been drawn on the basis of these pairs. It appears 
30 from Figures 7A and 78 that the selection and grouping according to the inven- 
tion transform the vague set of points shown in Figure 7A into two clear char- 
acteristic curves, of which one represents upward control (in the positive direc- 
tion) and the other downward control (in the negative direction). This simple 
example shows that the selection and grouping have great significance for the 
35 accuracy of the method according to the invention. 
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Calculation of characteristic curves 

[0034] In an embodiment after the selection and division, character- 
istic curves are calculated from the sample pairs in step 54. There are two 
characteristic curves, both of which are calculated by the same algorithm. The 

5 following describes an algorithm for calculating one characteristic curve. 

[0035] A control area (value range of control signal) is divided into 
bins uo(1)...uo(nbin), in which the values of the characteristic curve are stored 
as the calculation proceeds. Points Uo(1)...Uo(nbin) denote bin locations on the 
u-axis representing the control input and bin denotes the number of bins. Pa- 

1 0 rameters yo(1 ). . .yo(nbin) represent the values of output (measurement) y on the 
y-axis that represents the output. Thus the pairs uo(1)/yo(1)..-.Uo(nbin)/yo(nbin) 
formed by the bin locations and values define a characteristic curve. In addi- 
tion, the number of sample pairs used in the calculation of each bin value is 
calculated in counters nct(1)...nct(nbin). The number of bins may be any de- 

1 5 sired one. for example 3, 5 or 1 0. The bin Uo(1 ). . .Uo(nbin) values yo(1 )• • .yo(nbin) 
are updated by means of the selected sample pairs using a predetermined 

weighting function. 

[0036] In an embodiment of the invention, two pairs bp, bn-i are up- 
dated according to each sample pair and these pairs are selected so that the 
20 following condition is fulfilled: 

b„< ""'"""-^ (n,„-l)-H<6...; 

U —U- 
**inax nun 

where 

25 Um is the minute mean value of the control sample pair j, 

\/m is the minute mean value of the measurement of sample pair j, 
n = 1...bin, 

Umin and Umax are the minimum and the maximum of the control area, 
respectively. 

30 [0037] After this, the values yo(bn) and yo(bn+i) of the selected bins 

bn, and bn-i are updated as follows 
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where 

nct(1)...nct(nbin) each represents the number of updates (con 
trol/measurement pairs used) of each bin, 
Wn and Wn+1 are weighting coefficients 



W =1-J— 



l"m-"min-(^«-lK/| 



1 0 Ust is the distance between bins 



[0038] Finally, the numbers nct(bn) and nct(bn-1) of updates of the 
1 5 selected bins bn and bn-1 are updated as follows 
nct(bn) = a*nct(bn) + wi 
nct(bn+i) = a*nct(bn+i) + wa 
where a is a constant. 

[0039] Constant a is a 'forgetting factor'. Usually a = 1 is selected as 
20 the value of forgetting factor a, which means that forgetting is not in use. How- 
ever, in continuous monitoring it may be advantageous to select a value a < 1 . 
typically a = 0.9-0.9999. This value means that a newer measurement obtains 
a higher weighting coefficient in bin updating than an older measurement. Thus 
the influence of older measurements disappears gradually, i.e. they are "for- 
25 got". 

[0040] It should be noted that the method of calculating curves de- 
scribed here is only one feasible calculation method; yet it improves the accu- 
racy of calculation. Based on the above description, it is obvious to a person 
skilled in the art that other methods can also be used for calculating character- 
30 istic curves on the basis of the minute mean values selected and grouped ac- 
cording to the invention. 
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Hysteresis calculation 

[0041] The calculation 54 of characteristic curve produces two 
characteristic curves (or points of characteristic curves). One of the curves 
represents ascending control and the other one descending control. Hysteresis 
5 is calculated as a distance between the characteristic curves at one or more 
points in the direction of the control axis (u-axis). The selection of the number 
and location of the points preferably depends on how control has moved dur- 
ing a collection period. 

[0042] Hysteresis calculation will be exemplified at one point with 
10 reference to Figure 8. In the example, two points are known in one characteris- 
tic cun/e (ua, yA) and (uc, yc) and one point in the other characteristic curve 
(ub, ye) so that yA > ye > yc The distance between the characteristic curves in 
direction h of the u-axis is 

15 

h=u^-UB+ 

iyc-yA) 

[0043] Figure 9 illustrates two characteristic curves where hystere- 
sis has been calculated at several points. The values of the control axis (u- 
axis) of the characteristic curves shown in Figure 9 describe the control's rela- 
20 tive share (%) of the whole control area (100%). 

Example 

[0044] The calculation of the characteristic curve in the embodiment 
described above will be described here by an example. For the sake of clarity, 
the example describes only the calculation of one characteristic curve. In prac- 
25 tice, the hysteresis calculation updates two characteristic curves according to 
the same principle. 

[0045] Parameters used in the example 

Number of bins: nwn = 4 
30 Control min: Umin = 0 

Control max: Umax = 30 

Thus the distance between the bins is Ust = 10. 
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[0046] Before updating, the vectors describing the characteristic 
curve do not actually include any information. The vectors are 

• Measurements yO of characteristic curve 
5 • Controls uO of characteristic curve 

• Frequency information (number of hits) net on characteristic curve 

yO = [0 0 0 0] 
uO = [0 10 20 30] 
10 net =[0 0 0 0] 

[0047] The following pair of control and measurement is added to 
the characteristic curve 

15 Um = 12 

ym = 3 

First, bins are selected according to the following fomiula 

20 b,< K,„-i)-t-i<^>2 



(n^ - 1) + 1 = ^^-^ (4 - 1) + 1 = 2.2 => fti = 2, £»2 = 3 
Weighting coefficients are calculated 

10 

[0048] Measurement vectors yO are updated according to the fol- 



25 



lowing formula 

nct(b,)yo(b,) + w,y^ _ 0 + 08x3 _ ^ 
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[0049] Measurement vectors yO, bin 3 and bins 2 and 3 of the con- 
trol vector uO are updated in the same manner. 

Vector net is updated 

5 

nct(2) = a*nct(2) + Wi = 1*0 + 0.8 = 0.8 

and the same is perfomied on bin 3. Forgetting factor a=1, which 
means that forgetting is not in use. 

10 

[00501 After the first pair of control and measurement, the vectors of 
the characteristic curve are: 

yO = [0 3 3 0] 
15 uO = [0 12 12 30] 

net =[0 0.8 0.2 0] 

[0051] Next, the characteristic curve is updated using the following 
pair of control and measurement 

20 

Um = 24 

ym = 5 

[0052] Now the control is between bins 3 and 4 (bi = 3, ba = 4) so 
25 that it is closer to bin 3 (weighting coefficient Wi = 0.6) than to bin 4 (weighting 
coefficient W2 = 0.4). 

The measurement vector is updated as follows 

30 ^o^^^) nctib,) + w, 0.2 + 0.6 

The same fomnula is used in updating the control 

0.2x12 + 0.6x24 
"°^^>= 0.2 + 0.6 

35 
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vector net is updated, 

nct(3) = nct(3) + w1 = 0.2 + 0.6 = 0.8 

5 y0(4), u0(4), and nct(4) are updated as described above (the bin is 

empty before updating). Weiglitlng coefficient w2 = 0.4. 

[0053] After the second point, the vectors for the characteristic 



curve are 



10 



yO = [0 3 4.5 51 
uO = [0 12 21 24] 
net = [0 0.8 0.8 0.4] 

1 5 The following pair is still added 

Um = 15 

ym = 4 

20 The result is 

yO = [0 3.3846 4.3077 5.0000] 
uO = to 13.1538 18.6923 24.0000] 
net = [0 1.3000 1.3000 0.4000] 



25 



[0054] Assume further that the calculation has been continued by 
adding 97 pairs of measurement and control such that the control is a ran- 
domly selected integer between 0 to 30 and the measurement is the square 
root of this rounded to the nearest integer. 
30 [0055] The upper image in Figure 10 illustrates the pairs (x) of con- 

trol and measurement used in the example and the calculated characteristic 
curve (— ) yO vs. uO. The lower image illustrates the frequency curve net vs. 
uO (--^). 

[0056] In all the embodiments described above and in other vana- 
35 tions of the invention, the measurement y can be replaced with the set point r 
or a signal representing it. 
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[0057] The description and the related figures are only intended to 
illustrate the principles of the present invention. Various embodiments, varia- 
tions and modifications are obvious to a person skilled in the art on the basis of 
this description. The present invention is not limited by the examples but it may 
5 be modified within the scope and spirit of the attached claims. 



wo 2004/025384 — ^CT/FI2003/000669 



18 

CLAIMS 

1. A method of determining hysteresis of a process device in a 
process, the method comprising 

collecting sample pairs of signals representing a control input and a 
5 process output of a process device, 

calculating the hysteresis of the process device on the basis of the 

collected sample pairs, 

characterized by 

each collected sample pair comprising an average control input and 
10 an average output for a predetermined collection period , 

screening (53) the sample pairs suitable for further processing at 
least according to the magnitude and direction of the relative change of the 
average control input of each sample pair, 

grouping (53) the sample pairs into a first and a second group at 
15 least according to the direction of the relative change of the average control 
input of each sample pair, 

calculating (54) a first and a second control/output characteristic 
curve by means of the screened sample pairs of the first and the second 
group, respectively, 

20 detemiining the hysteresis of the process device as a distance be- 

tween the first and the second characteristic curve at least at one point. 

2. A method according to claim 1, characterized in that the 
signal representing the process output comprises one of the following: a meas- 
ured process variable (y). or a set point (r) of a process variable in process 

25 control. 

3. A method according to claim 1 or 2, characterized in that 

the screening (53) comprises 

calculating the change of average control input for each sample pair 
in relation to the average control input of the previous sample pair and the di- 

30 rection of the change, 

if the calculated change is smaller and has occun-ed in the same di- 
rection as the change calculated for the control input of the previous sample 
pair, the sample pair is selected for further processing; othenwise the sample 
pair is rejected. 
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4. A method according to claim 1,2 or 3. characterized in 
that the direction of the change is determined by calculating the sign of the 
change, in which case the sign is positive if the control input is increasing, and 
negative if the control input is decreasing. 

5. A method according to claim 1 , 2, 3 or 4, c h a r a c t e r I z e d in 

that the screening comprises 

selecting a sample pair UmG). YmO) for the calculation of the charac- 
teristic curve if the following two conditions are fulfilled: 

lAM„a)|<lA«.a-i)l 

sign{Lu^ (J)) = sigiiiLu^ (j - 1)) 



group, 



where 

UmG) is the average control input of sample pair j, 
ym(j) is the average measured output of sample pair j, 
15 j is an integer index, 

AUmG)= UmC)-UmG-1) 

sign(AUm) calculates the sign of Aum- 

6. A method according to claim 5, characterized in that the 
sample pairs are grouped into a first and a second group on the basis of the 

20 following rule: 

if the sign of Aum is positive, the sample pair belongs to the first 

if the sign of AUm is negative, the sample pair belongs to the second 
group. 

25 7. A method according to claim 1 , 2, 3 or 4, c h a r a c t e r i z e d in 

that the collecting comprises 

taking a momentary sample pair (ys, Us) from the control input (Um) 
and measured output (ym) at certain intervals, the interval being preferably in 
the order of one or more seconds, 

30 calculating the average sample pair (Um, ym) from the momentary 

sample pairs (ys. Us) taken during the collection period, which is preferably in 
the order of one or more minutes. 

8. A method according to any one of the preceding claims, char- 
acterized by calculating a reliability value for hysteresis by a function, 

35 which includes the following information as parameters 
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the number Ny1 of average sample pairs (Um, Ym) belonging to the 
first, ascending characteristic curve, 

the number Ny2 of average sample pairs (Um. Ym) belonging to the 
second, descending characteristic curve, 
5 the number Noy1 of sample pairs which belong to the first charac- 

teristic curve but are below the second characteristic curve, and 

the number Noy2 of sample pairs which belong to the second char- 
acteristic curve but are above the first control curve. 

9. A method according to claim 8, characterized in that 

^ 0 the reliability value L is calculated by function L=max(0.1 -Noy1 /Ny1 - 

Noy2/Ny2), in which case L may obtain values from 0 to 1, where L=1 is com- 
pletely reliable and L=0 is completely unreliable. 

10. A method according to any one of the preceding claims, 
characterized in that the calculation of each characteristic curve com- 

15 prises 

dividing the control area of the control input into bins uo(1)...uo(nbin). 
where uo(1)...uo(nbin) is the location of the bins on the u-axis representing the 
control input and bin is the number of bins, in which case the values 
yo(1)...yo(nbin) included in the bins represent output values on the y-axis de- 
20 scribing the output, and value pairs uo(1 )/yo(1 ) uo(nbin)/yo(nbin) define a char- 
acteristic curve, 

updating the values yo(1)...yo(nbin) of the bins uo(1)...uo(nbin) by 
means of the average sample pairs using a predetermined weighting function. 
1 1 . A method according to claim 10, characterized by 
25 updating two bins bn, bn-i according to each sample pair, the bins 

being selected so that the following condition is fulfilled 

^max ^min 

where 

Um is the minute mean value of the control of sample pair j, 
30 y^ is the minute mean value of the measurement of sample pair j, 

n = 1...bin, 

Umin and Umax are the minimum and the maximum of the control area, 
respectively, 

updating the values yo(bn) and yo(bn+i) of the selected bins bn, and 
35 bn.i as follows 
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10 



"-^ nct{b„) + w, 

where 

nct(1)...nct(nbin) each represents the number of updates (con 
trol/measurement pairs used) of each bin, 
Wn and Wn+1 are weighting coefficients 

w„ =1 



Ust is the distance between the bins 



max "min 



15 

updating the numbers nct(bn) and nct(bn-l) of updates of the se- 
lected bins bn and bn-1 as follows 
nct(bn) = a*nct(bn) + wi 
nct(bn+i) = a*nct(bn+i) + wa 
20 where a is constant. 

12. A system for determining hysteresis of a process device in a 

process environment, the system comprising 

means for collecting sample pairs of signals representing a control 
input and a process output of the process device, 
25 c h a r a c t e r i z e d in that each collected sample pair comprises 

an average control input and an average control output for a predetermined 
collection period, and the system comprises 

means (53) for screening sample pairs that are suitable for further 
processing at least according to the magnitude and direction of the relative 
30 change of the average control input of each sample pair, 
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means (53) for grouping the sample pairs into a first and a second 
group at least according to the direction of the relative change of the average 
control input of each sample pair, 

means (54) for calculating a first and a second control/output char- 
5 acteristic curve by means of the first and the second screened sample pair (Um, 
Vm), respectively. 

means (55) for determining the hysteresis of the process device as 
a distance between the first and the second characteristic curve at least at one 
point. 

10 13. A system according to claim 12, characterized in that 

the signal representing the process output comprises at least one of the follow- 
ing: a measured process variable (y), or a set point of a process variable (r) in 
process control. 

14. A system according to claim 12 or 13, characterized in 
1 5 that the screening means comprise 

means (53) for calculating the change of the average control input 
for each sample pair in relation to the average control input of the previous 
sample pair and the direction of the change, and for selecting a sample pair for 
further processing if the calculated change is smaller than the change calcu- 
20 lated for the control Input of the previous sample pair and has occunred in the 
same direction; othenwise the sample pair is rejected. 

15. A system according to any one of claims 12 to 24, c h a r a c - 
t e r i z e d in that the system comprises means (55) for calculating a reliability 
value for hysteresis by a function which includes the following infomnation as 

25 parameters 

the number Nyl of average sample pairs (Um, ym) belonging to the 
first, ascending characteristic curve, 

the number Ny2 of average sample pairs (Um. ym) belonging to the 
second, descending characteristic curve, 
30 the number Noyi of sample pairs which belong to the first charac- 

teristic curve but are below the second characteristic curve, and 

the number Noy2 of sample pairs which belong to the second char- 
acteristic curve but are above the first control curve. 

1 6. A system according to claim 14, characterized In that 
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the reliability value L is calculated by function L=max(0.1-Noy1/Ny1- 
Noy2/Ny2). in which case L may obtain values from 0 to 1 , where L=1 is com- 
pletely reliable and L=0 completely unreliable. 

17. A program product which includes a program code, which im- 
5 plements the method steps of a method according to any one of claims 1 to 1 0 
when it is run on a computer or the like. 
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